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It is believed that the presence of anticrossings with exponentially small gaps between the lowest 
two energy levels of the system Hamiltonian, can render adiabatic quantum optimization inefficient. 
Here, we present a simple adiabatic quantum algorithm designed to eliminate exponentially small 
gaps caused by anticrossings between eigenstates that correspond with the local and global minima of 
the problem Hamiltonian. In each iteration of the algorithm, information is gathered about the local 
minima that are reached after passing the anticrossing non-adiabatically. This information is then 
used to penalize pathways to the corresponding local minima, by adjusting the initial Hamiltonian. 
This is repeated for multiple clusters of local minima as needed. We generate 64-qubit random 
instances of the maximum independent set problem, skewed to be extremely hard, with between 
10^ and 10® highly-degenerate local minima. Using quantum Monte Carlo simulations, it is found 
that the algorithm can trivially solve all the instances in ~ 10 iterations. 



I. INTRODUCTION 

Adiabatic quantum computation (AQC) [T] is an im- 
portant paradigm for universal quantum computation 
[21 |3]. In a simple quantum adiabatic algorithm, the 
Hamiltonian of the system is written as 

H = A{s)Hb+B{s)Hp, 0<s<1, (1) 

with A{Q) > B{0) and A{1) < B{1). Here, s = t/tf is 
the dimensionless time with tf being the total evolution 
time. The system starts from a known ground state of 
Hb at i = and ideally ends in the ground state of Hp at 
t ^ tf, which is the solution to a problem. To ensure that 
the system ends up in the final ground state with high 
fidelity, the evolution should be very slow (adiabatic). In 
a closed system, the total computation time is related to 
the minimum energy gap gm between the ground state 
and first excited state of the Hamiltonian. For problems 
with very small gap at s = s* , a two-state approximation 
near the anticrossing yields the success probability Pf = 
1 — e"*^/'", with the adiabatic time scale |4j 

_ Ah {Q\dH / ds\l) 



where, |0) and |1) denote the ground and first excited 
states. To determine the computational complexity of 
AQC, therefore, one needs to know how gra scales with 
the problem size. This is indeed not an easy task, al- 
though analytic calculations have been possible for a few 
special examples [SHZ] . 

An important subset of quantum adiabatic algorithms 
is adiabatic quantum optimization (AQO), for which the 
final Hamiltonian Hp is diagonal in the computation ba- 
sis. The desired final ground state of the system is there- 
fore the classical global minimum of Hp. In this paper, 
we restrict ourselves to Hp and Hp of the form 

Hp = Y^ h^at + ^ J,,af a| , Hb = -Y. A.af , (3) 

i i,j i 



where erf and erf are Pauli matrices, and hi, Jij, and 
are real- valued parameters. 

One of the mechanisms that can result in a very small 
gm in AQO is when an eigenstate corresponding to some 
local minima of Hp (anti-) crosses that of the global 
minimum near the end of the evolution [8 . In this 
case, grn decreases exponentially with the Hamming dis- 
tance between the two minima. Perturbation expansion 
is proven to be a useful tool to examine these prob- 
lems [5HTT] , as these anticrossings typically happen when 
A = A{s)/B{s) ^ 1. To low orders of perturbation in 
A, the perturbed energy levels cross close to the point 
where an anticrossing exits in the exact spectrum of the 
system. For this reason, we call this type of anticrossings, 
perturbative crossings. They are also known as first or- 
der quantum phase transition points because of the non- 
analyticity of the ground state at these points in the ther- 
modynamic limit. The position of the antisrossings can 
be approximately predicted using low order perturbation 
expansion, while estimating gm requires higher orders [S] . 

Considering random exact cover instances, it was 
shown 9J that the probability of having such perturba- 
tive crossings increases with problem size, which was 
later numerically confirmed |12j . Others also found 
evidence of first order quantum phase transitions, even 
where the perturbation expansion is expected to break 
down [T5HT5] . These findings implied that NP-complete 
problems may not be solved efficiently using AQO. In 
most investigations on the scaling of AQO, however, 
at least some of the following assumptions are quite 
commonly made: 

1. Hp has no free input parameters. 

2. Ai is uniform among all qubits. 

3. All minima are non-degenerate states. 

4. Any suboptimal solution is an undesirable output. 

5. The computation is run exactly once. 

6. The system is completely isolated. 

Assumption 1 is not always the case, as several differ- 
ent Hp can represent the same problem. For Hp repre- 
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senting the Maximum Independent Set (MIS) problem, 
it was shown that a significant increase in gap size can be 
obtained by changing free parameters [TTl [TB] . Assump- 
tion 2 was examined in Ref. UF and it was found that 
choosing values randomly can have a non-zero proba- 
bility of eliminating the anticrossing between global and 
local minima. Assumption 3 was questioned by Ref. 1171 
who also criticized Ref. |9i based on neglecting the corre- 
lations between the local minima. Assumption 4 is not 
valid, since it is known that some optimization problems 
are NP-hard to solve approximately (TB]. It is also com- 
monly assumed that, at least for a closed system, there is 
no benefit of running the algorithm multiple times. For 
the case of open quantum systems, it has been suggested 
that running the evolution faster but multiple times may 
lead to much higher probability of success than running 
the evolution once but very slowly [H [19] . 

In a recent publication [llj, we examined the first 
5 assumptions for solving MIS problems via AQO. We 
showed that by changing the free parameters of the 
Hamiltonian, there always exists an adiabatic path along 
which no perturbative crossings would occur. However, 
no constructive way to find one of those paths was sug- 
gested. In this paper, we design a heuristic algorithm 
to determine such a path for a general problem. The 
objective is not to create a fiawless algorithm, but to 
demonstrate that a very simple heuristic algorithm, sim- 
ilar to adaptive simulated annealing [20 , can successfully 
eliminate severe anticrossings in hard instances. 



II. MAXIMUM INDEPENDENT SET PROBLEM 

Consider a general graph Q = {V,£), with a set of 
nodes (vertices) V and a set of edges £. An indepen- 
dent set in 5 is a set M C V such that no two of the 
nodes in M are adjacent. Every subset of an indepen- 
dent set is also an independent set. Naturally, there are 
independent sets that are not subsets of a larger indepen- 
dent set. We call them maximal independent sets. Thus, 
an independent set is maximal if all nodes outside the 
set are adjacent to at least one member of the set. The 
largest possible maximal independent set, which is also 
the largest possible independent set, is a maximum inde- 
pendent set. The problem of finding an MIS in a general 
graph is called the maximum independent set problem, 
and it is known to be NP-hard |21| . 

Suppose we associate a variable Xi € {0, 1} to each 
node. For a state x — [xi, X2, we define a cost 
function: 



£;(x) 



+ c 



E 



X.jX'j 



(4) 



with c > 1. Consider an independent set M of size to. 
Define a state x such that Xi = 1 if z S M, and Xi — 
otherwise. Since the nodes with Xi — 1 are not adjacent 
to each other, they don't contribute to the coupling c- 
term. In other words, since for every adjacent pair of 



nodes, at least one of the two has Xi = 0, the contribution 
of the coupling term to the cost function is 0. Therefore, 
the cost function will be E{x.) = Em = ~m. If M is a 
maximal independent set, it is easy to see that Em is a 
local minimum of £'(x). This is because, by removing a 
node from M, we will decrease ni and therefore increase 
i?(x), and by adding another node to it, we will increase 
i?(x) by at least c — 1 > 0. The latter is because the 
added node will be adjacent to at least one of the nodes 
inside M, due to the definition of maximal independent 
set. Therefore, maximal independent sets are the local 
minima of i?(x) and maximum independent sets are the 
global minima of £'(x). 

One can represent £'(x) in terms of an n-qubit Hamil- 
tonian by substituting x^ ^ {\ + <jf) /2. Ignoring a con- 
stant energy shift, we get Hp as in Q with: 



hi = - 



UiC ^ 



Jij — 



c/4 if ije£ 
otherwise 



(5) 



where Ui is the number of edges connected to (or degree 
of) the node i. 



III. PERTURBATIVE CROSSINGS 

For perturbation expansion, it is easier to work with 
the re-scaled Hamiltonian H = H/B{s) = Hp + \Hb, so 
that A = oo for s — and A = for s = 1. The eigen- 
functions of H and H are the same, but their eigenvalues 
differ by a factor of B{s). At the end of the evolution 
(A = 0), the eigenstates of H are the same as those of 
Hp. The ground state is therefore the global minimum 
of Hp and the low lying excited states are either local 
minima of Hp or states in the neighborhood of global 
or local minima. At small A, one can use perturbation 
expansion to calculate the eigenvalues and eigenstates of 
H. Perturbation expansion is valid as long as A < Ac, 
where Ac is the convergence radius of the expansion. 

Suppose there are K maximal (or maximum) indepen- 
dent sets Mk of size to, with k ^ 1, ...,K. States \Mk) 
and also every superposition of them are therefore degen- 
erate eigenstates of Hp, with energy E^^^ . Perturbation 
in A removes this degeneracy. Let \M)=J2k C^k\Mk) rep- 
resent the lowest energy superposition immediately af- 
ter the degeneracy is lifted. With the positive sign of 
Ai, all C'k will be positive real numbers with the con- 
straint: J2k — 1- Coefficients Ck can be obtained by 
partial diagonalization of H in the subspace of the local 
minima. The perturbed eigenvalue of this state can be 



written as: £'m(A) = E^^ + XE 



E 



i'^ = {M\Hb\M), and 
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(0) 



where 



(6) 



The perturbation Hamiltonian Hb causes single qubit 



fiips. The first order correction E 



(1) 

M 



0, because all 
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Mk are the same size (m), hence one cannot get from 
one minimum to another by a single bit flip (adding or 
removing a single node). The second order correction is 



E 



(2) 
M 



E 



k,i 



(7) 



where B^ i is the cost of flipping qubit i from state 
and the prime sign on the sum means that the sum is 
over all paths from \Mk) to |Mfe/) with two bit flips by 
flrst flipping qubit i and then qubit j . This also includes 
k — k' , which means flipping qubit i two times. Notice 

(2) • 

that for positive, real Ck and Ck', Ej^/ is always nega- 
tive for minima. Therefore, the second order perturba- 
tion correction always reduces the energy of eigenstates 
representing minima. 

Now suppose that M is the unique MIS of graph Q 
and represents the corresponding ground state oi Hp 

with eigenvalue e''^\ Also suppose there exist K max- 
imal independent sets M^, k = 1,...,K, with the same 
size m', producing degenerate local minima |Af^) of Hp 

with eigenvalue E^^]. Equation therefore, gives the 
perturbed energy of the above states (M — >■ M' for the 
local minima). The two states cross at A = A*, where the 
perturbed energies are equal: EmW = Em'W- Up to 
the second order perturbation, we have E^^j + A*^£'^^'' = 



E 



(0) 
M' 



A*^^^^^!, which leads to 



y-X^-iE7,-E-i>)/iE^,-E)^). 



(8) 



Since E^j^j, > E^^j\ in order for A* to have a real value 
we need E^^l < E^^ (or \E^^l\ > \E^^^\ since both cur- 
vatures are negative). This means that the local min- 
ima should have more negative curvature than the global 
minimum. The magnitude of the curvature in ([7| be- 
comes large if the energy cost B^^ i of bit flips from the 
local minima \Mk) is small. Moreover, if there are many 
degenerate local minima with two bit flip Hamming dis- 
tance from each other, each pair of those adds 4 terms to 
0. This means that if there is a large number of local 
minima connected to each other by 2-bit-flip paths, they 
may cause a large negative curvature creating an anti- 
crossing with the global minimum state. We call such set 
of nearby (in Hamming distance) local minima, a cluster. 
In practice, there could be many clusters of local minima 
and therefore there could be many anticrossings in the 
adiabatic path. 

Since the curvatures depend on according to ([T]), 
it could be possible to choose A^ in such a way that the 
inequality E'fp > E^fl^ would not be satisfied. Thus, 
the two states would not cross, at least up to second 
order perturbation. Our goal here is to construct a simple 
iterative algorithm that finds such A^ values heuristically. 



IV. AN ITERATIVE QUANTUM ALGORITHM 

Although it is difficult to determine nontrivial prop- 
erties of the global minimum, it is easy to gather infor- 
mation about the local minima. When there is a per- 
turbative crossing, if the quantum computation is run 
quickly (non-adiabatically) , with tf <C ta, the system 
wiU go to the excited state (|M')= X^fc C'fcl^fc))' after 
the anticrossing, instead of the the ground state {\M)). 
Repeating this process would sample the local minima 
|Mfc) with probability C|. 

We would like to use the above information to penal- 
ize the path to the corresponding local minima using 
Eq. (FtI). More specifically, we would like to reduce the 

— (2") 

curvature \E\p\ of state \M'), as much as possible, by 
changing A;, without significantly reducing that for the 
global minimum (|i?]j,j |). We cannot reduce all A^ to- 
gether, because it would simply reduce both. The bit 
flip energy costs Bk^i can be easily calculated once the 
local minima are known via the above sampling 

procedure. Determining all Ck, however, would require 
sampling all local minima, of which there could be expo- 
nentially many. Instead, we replace Ck' with Ck in ([7|, 
with minimal impact on the sum on average. We obtain 

- -E^» E |^a,.-5:a.^., (9) 



(fc,fc')j 



k,i 



Because is the probability of obtaining state |M^), we 
can compute fii by sampling results from running quickly: 



B 



k,i E Aj 



(10) 



sampled k 



From ([9]), it is evident that A^ with larger /x^ contribute 
more to the sum, and therefore should be reduced the 
most. What makes large is whether many of the de- 
generate local minima can be connected by first fiipping 
qubit i and then any other qubit j. When /^^ is small, 
Ai could be increased without significantly adding to the 
curvature, unless its effect on other /ij sums in ( 10 ) is 
large. 

For a given set of /i^, minimizing ^iMi? while keep- 
ing the geometric average dli^i)^^^ constant, yields 
Ai oc . When there is only one eigenstate \M') cross- 
ing the global minimum |A/), this choice of A.; might re- 
move the anticrossing even after the first iteration. How- 
ever, if there exists another state |M"), comprising an- 
other cluster of local minima |M^'), the above procedure 
may increase the curvature of Em", creating a new anti- 
crossing. This will take the system to a new set of local 
minima {\M^)) instead of the global minimum. The sec- 
ond iteration will penalize the path to without 
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TABLE I: The algorithm 



Initialize = 1 Vi. 

Anneal r times, saving each result. 

If a sufficient result ha s been returned, finish. 

Compute Hi using (10 1. 

Set Ai,ncw = A^^^^j^ , where /3 = + f) 
and K is the interation number. 
Rescale all A^'s to be within the feasible range. 
Go back to 2. 



considering \M'). This may now increase E^^p, which 
was reduced in the first iteration, and hence we will be 
back to the original set of local minima. Iterating such 
a process will only switch between those two clusters of 
local minima. In practice, there could be more than two 
clusters of local minima, so any successful algorithm must 
maintain a memory of the states reached in the previous 
iterations. 

This can be achieved by implicitly remembering pre- 
vious Ai values, in choosing A^^now oc Al~f^fi~^, i.e. a 
weighted geometric average with the previous value of 
A^. A smaller value of /3 will adjust less and remember 
more, so it may take more iterations to escape a partic- 
ular cluster of local minima. A large value of /?, on the 
other hand, makes the system prone to getting stuck in 
back-and-forth cycles between clusters of local minima. 
We found that f3 = where k is the current itera- 

tion number, will converge with high likelihood. Such a /3 
corresponds with performing a geometric average over all 
previous values of A^. Note that precise approximation 
of fii is not necessary to penalize the path to a cluster of 
minima, so a moderate number of samples are likely to 
be sufficient, and if not, the next iteration will build on 
this sample. 

The main algorithm is summarized in Table [T] In step 
5 of the algorithm, if some of Aj's are too small or too 
large, one can rescale them and limit A^'s that are still 
out of range to the acceptable minimum or maximum 
values. Note that step 4 can be completed in O(rn^) 
time on a single classical processor, or Q{\og{rn)) time 
with 0{rn^) classical processors, where r is the number 
of samples. 



V. TEST PROBLEM INSTANCES 

In order to adequately test this algorithm in simu- 
lation, we first required reasonably sized test problem 
instances with extremely small gm- Here, we focus on 
graphs with unique MIS (non-degenerate final ground 
state), as they represent harder instances than those 
with multiple MISs. Instances with a unique MIS are 
fairly uncommon in uniform random graphs. In testing 
80 graphs with 128 nodes and 1,572 uniform randomly 
placed edges, only 7 had a unique MIS, none of which 



had even remotely small minimum gaps, so generating 
uniform random graphs and hoping for instances with 
extremely small gaps would not have sufficed for gener- 
ating test instances. Generating random graphs from a 
distribution heavily skewed to have consistently small gm 
required a targeted approach. 

A key observation is that given a maximal independent 
set, randomly adding an edge between two of the nodes 
in the set makes the set dependent, and creates a degen- 
erate pair of maximal independent sets of 1 fewer node, 
separated by exactly 2 bit-flips. Doing this several times 
can produce many such pairs, which, as described in |11) . 
results in local minima eigenstates with large curvature 
as desired. 

Thus, the algorithm used for generating random 
graphs with extremely small gap anticrossings is: 

1. Create a graph with 64 nodes and 220 uniform ran- 
domly selected edges. 

2. Find, by a depth-first search through the space of 
independent sets, an independent set of size 20 to 
become the MIS, M. (The expected number of in- 
dependent sets of size 20 is approximately 5.7 mil- 
lion, so it is extremely likely that there is at least 
one.) 

3. For each node, i ^ M that is not adjacent to a node 
in M , uniform randomly select a node j £ M, and 
add an edge between i and j. This guarantees that 
M is a maximal independent set (not MIS yet). 

4. Continue the depth-first search until another inde- 
pendent set of size 20 is found. Remove one of the 
nodes and call this set M' . 

5. For each node i ^ A/' that is not adjacent to a node 
in M', if i ^ M, uniform randomly select a node 
j S M', and add an edge between i and j; if z S M, 
instead uniform randomly select j £ M' f] AI to 
ensure that M remains an independent set. 

6. Repeat 4 and 5 until no more independent sets of 
size 20 are found. 

The last step assures that M is an MIS. Caution should 
be used when applying this algorithm to much larger 
graphs, since although it can be executed in less than 
a second for 64 nodes, it does require time exponential 
in the number of nodes, assuming that the desired MIS 
size increases linearly with the number of nodes. 

All graphs generated with the above method are guar- 
anteed to have a unique MIS of size 20 (unless no inde- 
pendent set of size 20 was found in step 2). Of 51 graphs 
generated with this method, only 1 did not have an anti- 
crossing with gm small enough for the test scenario below 
(though it was still significantly smaller than that of all 
80 of the 128- node uniform random graphs examined). 
It was excluded from the test problem set, since it would 
be solved immediately. 
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FIG. 1: (Color online) A{s), the energy scale oi Hb, and B{s), 
the energy scale of Hp . 



All of these generated 64-node graphs have 100,000's 
of maximal independent sets (local minima), 1,000's of 
size 19 and 10,000's of size 18. Moreover, these tend to 
cluster into large groups of 1,000's of maximal indepen- 
dent sets connected by 2 bit-flip paths. Therefore, they 
are significantly harder than typical 64 qubit problems 
for AQO. 



VI. SIMULATIONS 

We would like to simulate this algorithm for the gener- 
ated problem instances. However, because for fixed-sized 
systems, "elimination" of a small gap anticrossing is not 
well-defined, in order to have a reasonable and objec- 
tive criterion for success, we must define time and energy 
scales. 

We choose as our A{s) and i?(s), energy scales ex- 
tracted from superconducting flux qubits similar in de- 
sign to those examined in [531 US] , which are plotted in 
Fig. [T] We choose r, the number of times the quan- 
tum computation is performed per iteration, to be 500, 
and the annealing time for each such computation to be 
tf — 0.08 ^s. Thus, if ta < 16 /xs, there is a high proba- 
bility (> 92%) that at least one of the 500 results from an 
iteration will be the global minimum. Also, for any result 
obtained that is not a minimum, one can easily perform 
gradient descent to reach a minimum, possibly the global 
minimum. Thus, if in an eigenstate crossing the ground 
state, the total probability of states that descend to the 
global minimum is > 0.005, there is a high probability 
(> 92%) that the global minimum can be obtained in 
this manner. Either case is considered to be successful. 

As described earlier, if the minimum gap is very small, 
such that tf <C ia, the system will occupy the excited 
state \M') after the anticrossing and one of the local min- 
ima will be reached, and the algorithm uses a sampling 
of these minima to calculate Hi using (10). Here, how- 




ever, we use Quantum Monte Carlo (QMC) simulations 
[T^ to provide the sampling needed. QMC provides sam- 
ples of computation basis states in the proportions that 



FIG. 2: (Color online) Visualization of perturbative crossings 
eliminated in 4 iterations. Each horizontal strip plots the 
ground state expectation {cri'^), for the ith qubit, as a function 
of s. There are 64 strips in each plot corresponding to the 
64 qubits. All qubits begin in a uniform superposition, so 
(ai'') = (green) at s = on the left. The final ground 
state represents the MIS, with -1-1 (red) for nodes in the set, 
and —1 (blue) for nodes not in the set. As s increases from 
to 1, the system localizes into low-energy minima (green 
moving toward red or blue). If it localizes into the global 
minimum, like the smooth transition in iteration 4, there is 
no perturbative crossing. However, if it localizes into local 
minima, there is at least one crossing, which will be visible 
as a sudden change in many qubits, as seen in iterations 0-3. 
Each iteration penalizes the path to the local minima into 
which the ground state localized, until no crossings remain. 
The penalization is strong enough that different local minima 
are found on each iteration. 



they appear in the ground state of the system. Before a 
small gap anticrossing, the ground state is approximately 
the superposition state \M') = J2k^k\Ml.). Therefore, 
using QMC to sample just before the anticrossing gives 
samples approximately as they would come from evolving 
the system. We use this fact to calculate Hi in (10 1. By 



incorporating gradient descent into the sampling, we can 
also determine the occupation of states in the well of the 
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FIG. 3: (Color online) Number of "unsolved" problem in- 
stances as a function of the number of iterations. All instances 
were solved within 13 iterations. 

global minimum, in addition to correcting for some of the 
single bit flip deviation from \M'). The adiabatic time 
([2]) can be computed using QMC in the same manner as 
described in [24j . If on an iteration, the adiabatic time 
is small enough to be computed accurately and is found 
to be < 16 /j,s, or the probability of states in the vifell 
of the global minimum just before a perturbative cross- 
ing is > 0.005, the instance is considered solved on that 
iteration. This is because the probability of finding the 
global minimum at least once in the 500 computations of 
the iteration is then > 92%. 

The "feasible range" of values, as mentioned in the 
algorithm description, was chosen to be between 1 /4 and 
8. More specifically, on each iteration, the values 
would always be scaled such that the smallest was 1/4, 
and the other A^ values rarely approached 8, especially 
after several iterations, where most A^ values were < 4. 

Beyond just testing whether the algorithm presented 
here is successful or not, it is critical to gain insight 
into how the changes in A^ values affect the evolution 
of the ground state. Although a 2^'^-dimensional sys- 
tem cannot feasibly be examined in detail, much can be 
seen by examining the ground state expectation values 

(a^z^) = (0|cri*^|0) of the 64 operators, i.e., the average 
magnetization of each qubit in the instantaneous ground 
state. Figure [2] illustrates these expectation values as a 
function of s for 4 iterations of the algorithm with an 
example instance. In each of the 4 plots, 64 horizontal 
strips color-code the value of (cri^^) for all 64 qubits, dur- 
ing the evolution, calculated using QMC sampling. At 
s = 0, the ground state is a uniform superposition of all 
computation states, so (cri*'') = (green), and at s = 1, 
the ground state is the global minimum of Hp, there- 
fore (ai^^) — ±1 (red/blue) depending on the state of the 
qubit. Moving away from s — 0, the expectations grad- 



ually tend toward —1 or -1-1 as the ground state settles 
into global or local minima. If the system directly settles 
into the global minimum (as in iteration 4), then there 

will be a continuous change of (cri''*) from to their ul- 
timate values with no sharp transition. If, on the other 
hand, the system initially settles into some cluster of lo- 
cal minima, then a sudden change is expected at the an- 
ticrossing between the eigenstates corresponding to the 
local and global minima (|M') and |Af)). Iterations 0-3 
clearly show such sudden transitions. After each itera- 
tion, the anticrossing moves to an earlier time, though 
this was not general among all the instances examined. 
In the last iteration, the anticrossing is completely re- 
moved and the transition from the beginning to the end 
state is smooth. 

All 50 of the extremely difficult 64-node random in- 
stances were solved within 13 iterations of the algorithm. 
The number of unsolved instances after each iteration is 
plotted in Fig. [3] Among the instances, 30 of them were 
solved in 2 iterations, and all but 1 were solved in 10 iter- 
ations, with an overall average of 3.0 iterations required. 
These data suggest that the selection of f3 to ensure equal 
application of penalties is quite robust at remembering 
previous penalties while still applying new penalties. 



VII. CONCLUSIONS 

We have demonstrated that a simple adiabatic quan- 
tum algorithm, based on penalization of paths to clus- 
ters of local minima by tuning singie-qubit tunnelling 
energies, is effective at eliminating extremely small gaps 
caused by perturbative crossings. We presented a method 
for generating 64-qubit random instances of maximum 
independent set with 10^ to 10^ highly degenerate local 
minima and a unique global minimum, causing pertur- 
bative crossings between the two. It is found that even 
for these instances, the algorithm can eliminate the per- 
turbative crossings in a small number of iterations. 
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